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The high frequency peaked BL Lac object PG 1553+113 underwent a flaring event in 2012. The High Energy 
Stereoscopic System (H.E.S.S.) observed this source for two consecutive nights at very high energies (VHE, 
E >100 GeV). The data show an increase of a factor of three of the flux with respect to archival measurements 
with the same instrument and hints of intra-night variability. The data set has been used to put constraints on 
possible Lorentz invariance violation (LIV), manifesting itself as an energy dependence of the velocity of light in 
vacuum, and to set limits on the energy scale at which Quantum Gravity effects causing LIV may arise. With 
a new method to combine H.E.S.S. and Fermi large area telescope data, the previously poorly known redshift 
of PG 1555+113 has been determined to be close to the value derived from optical measurements. 


1. Introduction 

PG 1553+113 is a high frequency peaked BL 
Lac object located in the Serpens Caput constel¬ 
lation. The object has been detected in VHE by 
H.E.S.S. m in 2006 and in high energies (HE, 
100 MeV<E<300 GeV) by Fermi [l|. The 7-ray 
spectrum presents the largest HE-VHE spectral break 
measured to date [H, @]. The source has an unknown 
redshift despite several attempts to measure it. The 
best estimate to-date, made by spectroscopy 1 , is 
0.43< ^ <0.58. 

In April 2012, the source underwent a flare reported 
by the MAGIC collaboration Q. Subsequently, the 
source has been observed by the H.E.S.S. telescopes. 
The data are used in this work to constrain its redshift 
and to probe a possible Lorentz invariance violation. 


2. Data Analysis 

2.1. H.E.S.S. data analysis 

The H.E.S.S. telescopes have observed 
PG 1553+113 in April 2012 for two nights. Data 
were analysed using the Model analysis Q with 
Loose cuts. The object has been detected with a 
significance of 21.5er in 3.5 hours of live time. The 
source spectrum is well fitted by a power-law model 
of the form: 

F{E) oc E+ 4 ' 85±0 ' 25 ) 

where E* = E /(327 GeV) and the flux is found to 
be 3.5 times higher than the measurements made in 
2005-2006 [|. 


Indications of intra-night variability have been 
found with the fit to a constant of the light-curve 
yielding a y 2 of 21.34 for 7 d.o.f. (P x 2 = 3.3 x 10 -3 ). 

Data taken in 2005-2006 were re-analysed and the 
spectrum is, in this case, well fitted by a log-parabola 
model: 

F(E) cx £-( 5 - 39±0 - 43 )-( 3 ' 95±1 - 40 ) l °SioE . 

where E* = E/(360 GeV). Both spectra (archival and 
flare) are presented on figure [0 

2.2. Fermi- LAT data analysis 

The Fermi- LAT data, from 300 MeV to 
300 GeV, have been analysed with the Science- 
Tools V9R32P5 and instrumental response functions 
P7REP_S0URCE_V15. A region of interest of 15 degrees 
has been used and the sky model has been built using 
the third Fermi catalog flO| . 

Data contemporaneous to the H.E.S.S. exposures 
taken in 2012 are well fitted by a power-law of index 
r = 1.72 ± 0.26. The pre-flare data are defined by the 
data taken from August 8, 2008 up to March 1st 2012. 
The measured spectrum is described by a log-parabola 
model (Fig. [IJ. 

Variability has been probed before, during and af¬ 
ter the flare using a bayesian blocks analysis )9]. No 
counterpart to the VHE flare was found in the HE 
light curve. 


3. Determination of the redshift 

The extragalactic background light (EBL) is a field 
of infrared photons that interacts with the VHE 7-rays 
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Figure 1: Spectral energy distribution in gamma-ray 
measured with Fermi- LAT and H.E.S.S. during the flare 
(red) and in the pre-flare state (see text) in blue. 


on their way to Earth. This absorption leaves a foot¬ 
print on the source spectrum that is used in this work 
to constrain its distance using a new bayesian model. 
The Bayes theorem reads P(0|F) oc P(0)P(V|0) 
where Y stands for the data and 0 the parame¬ 
ters. The likelihood, P(Y|0), is minimized during 
the spectrum determination. The model used here 
is a power-law corrected for EBL absorption i.e </> = 
N x (E/E 0 )~ r x e~ r( ' E,z \ and then 0 is IV, T and z. 

To construct the prior P(0), the following assump¬ 
tions have been made: 

• EBL-corrected power-law cannot be harder than 
the Fermi measurement, the prior being then a 
truncated Gaussian of mean T Fermi and width 
u that accounts for statistical and systematic 
uncertainties of both instruments. 

• Softening of this power-law, that arises from 
emission effects, is permitted with a constant 
probability. 

• It is also assumed that distant sources are harder 
to detect and that P(z) oc exp(—r(z)). 

The EBL absorption is computed using the model of 
[fj]. Marginalising over the parameters gives a redshift 
z = 0.49 ±0.04 (Fig. [2j) in agreement with other mea¬ 
surement 0] or constraints derived using GeV-TeV 
data j||. 


4. Lorentz invariance violation (LIV) 

Tests for a possible LIV effect were performed by 
searching for a non-zero dispersion parameter r ra (~ 
(AE ) n ) the H.E.S.S. data of the flare. This is done 
by testing a correlation between arrival times of the 
photons and their energies. 


A maximum likelihood analysis based on [?|, has 
been modified to tackle the non-negligible background 
present in the data. 

For noN events recorded in the ON-source region 
with arrival times ti and energies £), the likelihood 
reads: 

L(r n ) = niLT P{Ei,ti\T n ) with: 

P (Ei , U | r n ) = w s • P Sig (Ei , ti | r n ) + (1 - w s ) • P Bkg (Ei , ti) 

The probability Psi B was mainly determined from 
a parametrization of the light curve at low energies 
parametrization while P Bkg was built assuming a con¬ 
stant background. The factor w s accounts for the rela¬ 
tive weight of signal events with respect to background 
events. 

Constraints on r ra led to lower limits on the Quan¬ 
tum Gravity energy scale Eqq. The 95% 1-sided 
lower limits for the subluminal case are: Eqg,i > 
4.32 x 10 17 GeV and Eq G , 2 > 2.11 x 10 10 GeV for 
linear and quadratic LIV effects, respectively. Fig¬ 
ure 3 compares these results with other limits from 
less distant AGN flares. While the statistics is more 
limited here, the distance of the source makes the sen¬ 
sitivity to possible LIV effects comparable to previous 
results. 


5. Conclusions 

The VHE emitter PG 1553±113 underwent a flaring 
event in VHE with an increase of its flux by a factor 
of 3.5. No counterpart of this flare was found in the 
HE regime by Fermi. 

This data set has been used to constrain the red- 
shift of the source to be z = 0.49 ± 0.04 using a 
novel method based on 7 -ray data. The flare is 
also used to put lower limits on the LIV effect with 
Eqg.i > 4.32 x 10 17 GeV and E QG ,2 > 2.11 xlO 10 GeV 
for linear and quadratic effects. 
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Figure 2: Posterior probability for the redshift determination and comparaison with other measurements. 
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Figure 3: Lower limits on Eqg.i from linear dispersion (left) and on Eqg,2 from quadratic dispersion (right) for the 
subluminal case obtained with AGN as a function of redshift. 
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